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ABSTRACT 

We analyse the overionisation or proximity zone of the intergalactic matter around 
high-redshift quasars in a cosmological environment. In a box of 64 Mpc base 
length we employ high-resolution dark matter only simulations with 1024^ particles. 
For estimating the hydrogen tempe r ature and density distribution we use the effective 



equation of state bv lHui fc GnedinI (|1997[ ). Hydrogen is assumed to be in photoionisa- 



tion equilibrium with a model background flux which is fit to recent observations of the 
redshift dependence of the mean optical depth and the transmission flux statistics. At 
redshifts z = 3, 4, and 4.8, we select model quasar positions at the centre of the 20 most 
massive halos and 100 less massive halos identified in the simulation box. From each 
assumed quasar position we cast 100 random lines of sight for two box length includ- 
ing the changes in the ionisation fractions by the QSO flux field and derive mock Lya 
spectra. The proximity effect describes the dependence of the mean normalised optical 
depth ^ = Toff_ QSO /''cff, Lyct as a function of the ratio of the ionisation rate by the QSO 
and the background field, u) = FQso/ruvB, i-e. the profile C = (1 + w/a)~°-^, where a 
strength parameter a is introduced. The strength parameter measures the deviation 
from the theoretical background model and is used to quantify any influence of the 
environmental density field. We improve the statistical analysis of the profile fitting in 
employing a moving average to the profile. We reproduce an unbiased measurement of 
the proximity effect which is not affected by the host halo mass. The scatter between 
the different lines of sight and different quasar host positions increases with decreas- 
ing redshift, crioga ~ 0.08, 0.20 and 0.36 for z = 4.8, 4, and 3, respectively. Around the 
host halos, we find only a slight average overdensity in the proximity zone at comoving 
radii of 1 < rc < 10h~^ Mpc. However, a clear power-law correlation of the strength 
parameter with the average overdensity in Tc is found, showing an overestimation of 
the ionising background in overdense regions and an underestimation in underdense 
regions. 

Key words: diffuse radiation - galaxies: intergalactic medium - galaxies: quasars: 
absorption lines 



1 INTRODUCTION 

The Lyman alpha (Lya) forest in high-resolution quasar 
(QSO) spectra blueward of the QSO Lya emission repre- 
sents an excellent tracer of the high redshift matter distri- 
bution. According to det ailed sim ulations and a large obser- 
vational material (cp. e.g. lMeiksrni.2009. ) the forest of lines is 
generally ascribed to the absorption by a tiny fraction of re- 
maining neutral hydrogen H I. This H I forms the cosmic web 
at low gas overdensities of a factor of a few and arises natu- 
rally due to gravitational instabilities (|Petitiean et al ] | 19951 : 

* E-mail: apartl@aip.de (AMP) 



iMiralda-Escude et al. I ll996l : iHui et al.lll997l ) in a standard 
ACDM cosmological model. The intergalactic gas distribu- 
tion and its ionisation state forms from a balance of the 
ionising intergalactic UV radiation stemming from galax- 
ies and quasars and the developing inhomogeneous density 
field. Due to the high ionisation degree, the mean free path 
of UV photons in the Lya forest is hundreds of Mpc long. 
Hence, many sources which are distributed over large dis- 
tances contribute to the UV background (UVB), creating a 
quite homogeneous UVB. 

This picture changes slightly in the vicinity of a QSO. 
There the QSO radiation dominates over the overall back- 
ground, additionally increasing the ionisation state of the 
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intergalactic gas. This manifests itself in a reduction of 
the number of observed absorption features towards the 
emission redshift of the QSO. This proxi mity effect was 
first o bserved in low resolution spe ctra by Carswell et all 



1986) 



1983) 



l|l98dl_ and was later confirmed by Murdoch et alj 
and Irvtleil (|l987l l. In a seminal paper, Baitli k et al. 
used the proximity effect to estimate the intensity of the 
intergalactic UV radiation in comparison to the QSO lu- 
minosity, assuming that QSOs reside in random regions 
of the universe. Investigations of the proximity effect 
along these lines using absorption l i ne counting statistics 
dLu et all Il99ll IWihiger et al'l Il994l: ICristiani et"al] Il995l: 



Giallon go et al Ill996l: Isrianand fc Khardll996l : ICooke etall 
_1997; Scott et al.l |200ol'l or using pixel statistics of the 
transm it ted flux (ILiske fc W illigcr 200 ll : iDall'Aglio et all 
l2008bl R I2OO9I : ICalverlev et al.. JOllj) reveal the proximity 
effe ct in many high- r esolutio n QSO spectra. Recent results 
by IP airAgli o et all ll2008ah , ICalverlev et all (|201lh . and 
iHaard t & M adaul ( 201ll 'l indicate a steady drop in the ion- 
ising background radiation's photoionisation rate Fuvb 
lO"" '^ s"^ at 2 ~ 2 to TuvB ~ iQ-i^ ^s g-i for z > 5.5. 
However results obtained from SDSS spectra indicate a con- 
stant ionising ba ckground level between 2.5 < z < 4.6 
l|DaU'Aglio et al.| [2009). 

Independent measurements of the UV back- 
ground radiation can be obtained with the flux decre- 
ment method, which requires cosmological simula- 
tions to model sta ti stical propert i es of t he Lyg forest 
jRauch et all Il997l: iTheuns et all I l998l: I Songaila etlT 



19991: 



2003 



McDonald fc Miralda-EscudelTOOl: Mciksin fc Whit^ 



2005 



Tvtler et al1l2004 iBolton e t al. 200 5,: iKirkman et all 



Jena et all l20(35l : iFauc hcr- Gig uere et al.l I2OO8I : 



Bolton fc Haehneltl [20071 ). It becomes clear that 



measure- 



ments using the proximity effect seem to overestimate 
the UV background by a factor of a few. It is discussed 
whether this discrepancy arises from environmental effects 
such as clustering, gas infall, or the large scale density 
environmen t, which are neglecte d in the proximity effect 
modelling llRoflinde et al.l l2005l: ICuimarkes et al.i l2007l : 



iFaucher-Giguere et all 12003 : iPartl et al.ll2010l ). Such effects 
can result in an o veres timation of the UV background by 
up to a factor of 3 jLo eb fc Eis ensteinlll995l '). 

QSOs are thought to reside in very massive halos. Ob- 
servational estimates on the QSO halo mass revea led val- 
ues of a co uple of 10^^ (| d a Angela et~al] I2OO8I ') up to 
~ IO^Mq (|Rollinde et aljboOSi '). From numerical simula- 
tions of structure formations it is known that such massive 
halos are not located at random positions in the universe, 
but form in dense environments where galaxies cluster. Ob- 
servational determinations of the large scale environment 
around QSOs from Lya forest spectra indicate that QSOs 
are embedded in l arge scale ove rdensities extending from 
proper ^ 3-5 Mpc jP'Odorico e t al. 2008) to proper ~ 10- 
15 Mpc (RoUinde et al.' '2005'; 'Guimaracs ot al. 200?). Nu- 
merical simulations by Fauchcr-Gigucrc ct al. (2008) show 
large scale overdensities of proper ~ 3 — 6 Mpc for redshifts 
A > z > 2, consistent with the results by iD'Odorico et al.l 
l|200St ). Using detailed radiative transfer simulati ons of three 
QSOs residing in different cosmic environments [Partl et all 
hereafter PI) found indications that such large scale 
overdensities weaken the apparent proximity effect signal, 
resulting in an overestimation of the UV background. How- 



ever, the low number of QSO hosts in the PI study does 
not allow us to securely confirm the effect. We therefore ex- 
tend this study using a large sample of different dark mat- 
ter halos in various mass ranges and inquire how such large 
scale overdensities affect UV background measurements. It 
is also checked whether a possible dependence of the UV 
ba ckground with the host halo m ass exists, as was suggested 
bv lFaucher-Giguere et all (|2008l ). 

This paper is structured as follows. In Sect.[2]we present 
the dark-matter simulation used for this study and de- 
termine realistic models of the Lya forest from a semi- 
analytical model of the intergalactic medium. In Sect. [3] we 
introduce the proximity effect as a measure of the inter- 
galactic ionising background flux and characterise the halo 
sample used in this study. It is further discussed how Lya 
forest mock spectra are generated. Subsequently in Sect.|4] 
we evaluate the effects of large scale overdensities and infall 
velocities on the proximity effect. A possible dependence of 
the UV background measurements on the halo mass and on 
the large scale mean density around the QSO are assessed. 
We summarise our results in Sect.[Sl 



2 SIMULATION 

2.1 Distribution of baryons in the IGM 

In order to obtain a model for the gas content in the IGM 
that is in agreement with observed properties of the Lya 
forest, we employ the 6 4 Mp c[^ DM simulation of the 
CLUES projec10 (Gottl ober et al.|[201Q l with 1024=^ DM par- 
ticles. The GADGET2 (|Springel2005l ') simulation has a mass 
resolution of mp.pM = 1-86 x lO'^ Mq and uses a 

WMAP5 (|Hinshaw et al.ll2009l ) cosmology. 

The distribution of DM particles was obtained at four 
different redshifts z — 4.8, 4, 3, and 2. Using triangular 
shaped cloud (TSC) density assignment we obtain a 800'^ 
regularly spaced density and velocity grid from the DM par- 
ticle distribution. Halos have been identifled in the simula- 
tion u sing the hiera rchical friends-of-friends (HFOF) algo- 
rithm ([Klvpin et al.lil999 ) with a linking length of 0.17. 

To obtain an IGM gas density held from which we 
can construct Hi Lya forest spectra, it is assumed that 
the properties of the baryonic component, such as density 
and bul k velocity, are proportional to those of the dark 
matter l|Petitiean et al.lll995l : [Meiksin fc Wh"it3l200ll ). We 
have checked this assumption with a comparable GAD- 
GET2 gas dynamical simulation using 2 x 256'^ particles 
with a size of 12.5 h-^ Mpc (For ero-Romero et al 2011, 
m preparation). This corresponds to a mass resolution of 
'Tip,bar. ~ 3.7 X 10^ Mq in gas. The gas dynamical 
SPH simulation includes also radiative and compton cool- 
ing, star formation, and feedbacks throug h galactic winds 
using the model of ISpringel fc Hernguistl (i2003) . Further- 
more, a UVB generated from QS Os and AGNs is included 
at 2; < 6 l|Haardt fc Madaulll996l ). The state of the simula- 
tion has been recorded at 2 = 5 and 4. The simulation was 
rerun with DM only using the same realisation of the initial 



^ Distances are given as comoving distances unless otherwise 
stated. 

^ http://www.clues-project.org 
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conditions. We obtain density fields from tfie two simula- 
tions with equal spatial resolution as the larger 64 Mpc 
sized fields. The DM densities are then compared to the den- 
sities derived from the gas dynamical simulation and the re- 
sulting DM to gas density relation is shown in Fig. [T] The 
baryonic component follows the DM density closely, espe- 
cially at 2 = 5. For lower redshifts a linear relation is still 
obtained, however the scatter around the line of equality, 
where the DM overdensity is equal to the gas overdensity, 
increases. A small bump in the DM to gas density relation 
towards higher gas densities develops with decreasing red- 
shift at gas densities where a fraction of the gas becomes 
shock heated. The density in these shocks is slightly larger 
than the underlying DM, giving rise to the small bump in 
the DM to gas density relation. Studying the density proba- 
bility distribution reveals that a,t z — 5 the two components 
are similar (see Fig. [1]). At lower redshifts, the baryon dis- 
tribution is slightly shifted to higher densities and decreases 
steeper than the DM one at high densities. For completeness 
we additionally show the effective equation of state derived 
from the density and temperature fields in Fig. [T] 

From the gas dynamical simulation we find the assump- 
tion that baryons follow the DM distribution to be reason- 
able for the low density regions from which the Lya forest 
arises. We therefore use the DM only simulations to derive 
the gas density and velocity fields as described in PI. By 
using a TSC mass assignment scheme, we implicitly smooth 
our density field on sizes of 1.5 cells (i.e. 120 h'^ kpc). 
This is comparable to a constant Jeans length smoothing 
of Ajeans oc S'^^^il + z)'^!'^ ~ 150 K'^ kpc at 2 = 3, where 
5 is the density contrast. 



2.2 Model and calibration of the intergalactic 
medium 

The characteristics of the IGM gas are determined analo- 
gously to the method described in PI. We will therefore 
only briefly sketch the method with which a representative 
IGM model is obtained. 

The temperature in the IGM can be approximated 
by th e so called effective equation of state (|Hui fc GnedinI 
liggTl ). which is expressed as T = ro(l -|- J)'^"^ The effec- 
tive equation of state shows a linear trend for overdensities 
0.1 ^ 5 10 (compare Fig. [T|). For higher densities its be- 
haviour is approximated by assuming a constant value above 
a temperature cut-off Tc ut-og = T ( & = 1 0) . By employing 
the method developed in iHui et all (|l993), we compute H I 
Lya absorption spectra from the DM density and velocity 
fields, once the parameters To, 7, and the UV background 
photoionisation rate (Fuvb) are determined. These are con- 
strained by matching statistical properties of the simulated 
H I absorption spectra with observed relations. Such a cali- 
bration will be important in obtaining realistic representa- 
tions of the IGM that closely mimic the observed character- 
istics. 

In order to calibrate our IGM model, we em- 
ploy four observational constraints sorted by increas- 
ing importance: (i) The observed equation of st at e 
l|Ricotti et alJl2000l : ISchave et al.ll2000l : lLidz et al."20ld'). (h) 
the evolution of the UV background photoi o nisation rate 




1.5 2.0 2.5 3.0 



3.5 
z 



4.0 4.5 5.0 5.5 



dHaardt fc Madaul 12001 : Bianc Met akl I2OOII : iBolton et al] 
I2OO5I : iDaH'Aglio et al.ll2008h . i20ol)77iii) the observed evo- 



Figure 2. Upper panel: The evolution of the UV background 
photoion isation rate in the f our snapshots (black squares) com- 
pared to iBolton et al.l ||2005| ) (grey poi nts shifted by Az = 0.05 
for better visibility), meas urements by iDall' Aglio et al.l l l2008al ') 
(grey open diamonds) and ICalverlev e t al.l ll2 01ll) ( gr ey s tar), a 
fit to SDSS data measurements bj'jDair Aglio et al. (2009') (grey 
dashed Une) and predictions by ,Haardt &; Madau ( 2001) (grey 
line). Lower panel: The effective optical depth o f our models 
(black points) in comparison to measurements by ISchave et al.l 
(200^) (grey poi nts). The c o ntinu ous line shows the fit to obser- 
vational data bv lKim et al. I 1I2OO7I '). 



lutio n of the effective optic al depth in the Lya forest 
(Schave et al.ll2003l : iKim et al. 2007). a nd ( iv) the transmit- 
ted fiux probability distribution fFPD. lBecker et al.l[2007t ). 

The statistical quantities of the simulated spectra are 
derived from the simulation using 500 lines of sight ran- 
domly drawn through the cosmological box. Due to the avail- 
ability of high resolution spectra with high signal-to-noise 
(S/N) ratios of up to 120, we approximate such high quality 
data by assuming noi se free spectra. It has been noted by 
ICalverlev et al.l (|201ll ) that low S/N levels introduce a sys- 
tematic shift in the proximity effect signal. Since we want 
to quantify the physical effect of the halo's surrounding en- 
vironment on UVB measurements, adding noise would only 
lead to degeneracies between the two effects. However we 
will briefly address the influence of noise in Section [4.21 

The spectra are convolved with the instrument proflle of 
the UVES spectrograph and are then binned to the typical 
resolution of UVES spectra of 6.7 km s~^. We further as- 
sume that the QSO continuum can be perfectly determined 
and we therefore consider no uncertainties in the continuum. 
This procedure is adopted to precisely quantify physical ef- 
fects without contamination of the signal with observational 
uncertainties. These will only increase the variance in the re- 
sults discussed below. 

To determine the model parameters for the IGM, we 
adopt an iterative method based on minimisation be- 
tween simulated and observed FPD, as the latter provides 
one of the strongest constrains on the Lya forest properties. 
Our steps are as follows: 

(i) We choose a initial guess (Tq, 7, Fuvb) „ according to 
the measurements of the equation of state by ISchave et al.l 
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Figure 1. The top row shows the effective equation of state of a 12.5 Mpc gasdynamic simulation binned to cells of 80 kpc at 
redshift z = 5 (left panels) and z = 4 (right panels). The middle row shows the gas overdensity 1 + (5gas as a function of dark matter 
density 1 + 5dm- The white diagonal line marks the linear relation if gas strictly follows the dark matter. In both rows the colour coding 
refers to the logarithm of the respective number of cells. The lowest row shows the probability distribution of the dark matter (solid line) 
and gas (dashed line) overdensity. 
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Table 1. Model parameters of the semi-analytical model. 
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Figure 3. Upper panel: Comparison of our choice s for 7 (black 
square s) with observationally derived results by ISchave et al.l 
BOOOl ) (grey points). Lower panel: C omparison of o ur mo del Tq 
(black poin ts) with ob servations by ISchave et al.l l|200(t ) (grey 
points) and iLidz et al.l Jioig) ( grey open diamonds). 



ll2000l) andlLidz et al.l(l2010l ) and the evolution of the UVB 
bv lHaardt fc Madaul (|200ll ) (see Figs. [2] and [3)1 . 

(ii) With the above initial guess and using 500 lines 
of sight, we determine the average effective optical depth 
T"eff(z) = — In {F{z)) with F being the transmitted flux and 
the averaging is performed over the whole line of sight. 
As To and 7 have a subdominant effect on Tcs{z), we 
first tune Fuvb to obtain a match in the effective optical 
depths of the mo del spectra with recent observations from 
iKim et al.1 (|2007l ). This results in a new set of parameters 
(To,7,ruvB)i. 

(iii) Finally we construct our simulated FPD. As obser- 
va tional constra i nt we employ the log normal fits obtained 
bv lBecker et~aLl l|2007l ) in a redshift interval of Az — ±0.25 
centred on the snapshot redshift z (see Fig. |4]). By using 
their fits to the FPD, our model spectra can be compared 
with the observations without considering the effect of de- 
tector noise. However for consistency we also compare the 
si mulated FPD with th e raw combined observational results 
by HectereEai] l|2007l ). folding the simulation result with 
the global noise function of the combined observed sample 
(see Fig. [5]). The best fit parameters (To,7,ruvB), are then 
iteratively determined with a Simplex optimisation proce- 
dure. 

Our best fit parameters are presented in Table [TJ and 
plotted in Figs. [21 El SI and O in comparison with different 
literature results. 

The inferred evolution of th e UV background in Fi; 
[2] closely follows recent results bv lHaardt fc Madaul l|200 
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iBolton et all ()2005h . and lPaU'Agho et al.l (j2008b[ ). The val- 
ues for redshi ft z = 4.8 is rather high when compared 
with Calverl ev et al.l ()201ll ). however it is still in agreement 
within the 2a limits. Furthermore the effective optical depth, 
also show n in Fig. [2l is c onsistent w ith high resolution obser- 
vations bv lSchave et al.l (2003) and lKim et"all ()2007l ). Addi- 
tionally the final parameters of the equation o f state To and 
J lie w ithi n the measuremen t uncertainties of ISchave et al.l 
(|2000l) and lLidz etHI l|2010l) . which is evident from Fig.O 
The fiux probability distributions estimated from the 
simulated sight lines agree reasonably well with the noise- 
corrected observed profiles estimated bv lBecker et al.l (|2007l ) 
who assumed a log-normal distribution of the optical depth 
in the Lya forest (see Fig. [JJ . At redshifts z = 2 and 3 the 
match between the mean observed profiles and our simu- 
lations is good. However at higher redshifts, the agreement 
marginally decreases, even though the distributions are con- 
sistent within the variation between different lines of sight. 
This discrepancy manifests itself clearer when comparing 
the simulations with the raw observational data in Fig. [5] 
which include the effects of noise. The discrepancy is espe- 
cially strong at the high transmission end of the distribution. 
This is most certainly caused by our crude noise modelling 
of a sample with inhomogeneous S/N, which we obtained by 
stacking the noise functions of the various lines of sight. On 
the other hand we assume a perfect knowledge of the contin- 
uum level in the spectra, which is a challenge to determine 
in observed spectra, especially at high redshifts. 



3 METHOD 

3.1 Line-of-sight proximity effect 

Bright UV sources, such as QSOs, alter the ionisation state 
in their vicinity strongly up to proper distances of a couple 
of Mpc. The resulting change in the ionisation state directly 
translates into a decrease in optical depth. This decrease 
manifests itself in a reduction of the absorption line den- 
sity in the Lya forest when approaching t he redshift of the 
source. Using line-counting statistics, Bai tlik et al.l (|l988l ) 
measured the line-of-sight proximity effect in QSO spectra 
for the first time. They assumed that QSOs are situated 
in a mean IGM environment, neglecting any redshift dis- 
tortions from peculiar velocities of the absorption clouds. 
They further assumed that the QSO's radiation field is ra- 
dially decreasing from the source proportional to , i.e. 
geometric dilution. Using detailed radiative transfer simula- 
tions, PI showed that radiative transfer effects only play a 
marginal role in the line-of-sight proximity effect and that 
the assumption of geometric dilution holds. 

Assuming photoionisation equilibrium, the change in 
the optical depth as a function of distance r to the QSO 
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Figure 4. Probability distribution function of t ransmitted flux F at redshifts 2, 3, 4, and 4.8. Grey lines give the mean observed PDF 
obtained from the fits excluding noise given in ijBecker et al.ll2007t ) with 1(7 error bars denoting the variation between different lines of 
sight. The mean is derived in redshifts bins of ±0.25 centred on the model redshift. The black full line shows the results from our models. 



can be expressed as 

TPE(r) = rLy.M(l + ^M)"' (1) 

IILiske fc Willigei|[200ll ). where rpE is the observed optical 
depth, TLya is the optical depth in the absence of a QSO, 
and 



rQso(r) 



1 



(2) 



FuvB 167r2 Ji^ll.uvb 3 — Qq 

acts a s a normalised distance to the QSO ijPairAglio et al.] 
l2008bl ). Here rQso('') and Fuvb are the photoionisation 
rates of the QSO and the UV background respectively, 
L^LL is the QSO's luminosity at the Lyman limit, while 
J'l^LL.uvB is the UVB's Lyman limit flux, is the slope 
of the UVB's spectral energy distribution oc ly''^ , and Oq 
denotes the spectral slope of the QSO's emission oc u"'^. 
In this work we assume for simplicity the UVB to be 
domin ated by QSOs. We therefore follow [Saardt fc Madaul 
l|l996r) and chose the slope of the UV B to be Ob = —1.5 . 
This is consistent with observations by iTelfer et al.l (|2002| ) 
who find Qq — —1.57 ± 0.17. The effect of differing spectral 
shapes between t he UV B and the QSO has been studied in 
iDall'Aglio et all l|2008bl ') and PI. The proximity effect is in- 
troduced into lines of sight drawn from our simulation boxes 



by modifying the neutral hydrogen fraction in real space us- 
ing nHi,PE(r) = nHi,Lya(r) (1 -I- uj{r))~^. 

3.2 Environments 

To study the impact of the large scale environment on the 
proximity effect measurement as a function of QSO host halo 
mass, we pick halos in given mass ranges from the simulation 
to serve as QSO hosts. Around each halo, 100 lines of sight 
with a length of two comoving box sizes (128 Mpc) are 
randomly drawn from the box assuming periodic boundary 
conditions. Our sample of QSO host halos cover a wide range 
of masses. At each redshift we use the 20 most massive halos, 
50 halos with a mass around IO^^A^q, and 50 halos with 
a mass around 10^^ Mq. For the 10^^ and 10^^ Mq mass 
objects we choose 50 halos from the mass sorted halos with 
a mass larger than the cutoff. For the more massive bin, less 
than 50 halos with mass higher than the cutoff are present 
in the simulation at 2 ^ 4. In this case we choose all the 
halos above lO^'^ Mq and extend the range to lower masses 
until the sample consists of 50 halos. For these redshifts, 
the IO^^Mq mass range overlaps with the 20 most massive 
halos. The covered mass intervals are given in Table [2] as a 
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Figure 5. Pro bability distribution function of transmitted fiux F at rcdshifts 2, 3, 4, and 4.8. Grey lines with Icr error bars give the 
observed PDF ijBecker et al.l [20071 ) and the black full line show the results obtained from our models folded with the signal to noise 
function of the combined observations. The combined observational data is derived in redshifts bins of ±0.25 centred on the model 
redshift. 



Table 2. Mass range of the halo samples as a function of redshift. 
The 10^^ Mq and lO^^ Mq samples contain 50 halos each. The 
20 most massive halos are given in units of lO^^ Mq, the other 
mass bins in units corresponding to the mass bin. 



z 


20 most massive 


10^2 Mq 


10" Mq 


4.8 


2.07-0.92 


2.07-0.53 


1.04 - 1.00 


4 


4.84 - 1.51 


4.84 - 0.88 


1.03 - 1.00 


3 


15.8 ~ 3.41 


1.38 - 1.00 


1.01 - 1.00 


2 


40.2 ~ 6.19 


1.14 - 1.00 


1.01 - 1.00 



function of redshift. We further note that the halos in the 
10^^ Mq mass bin are certainly not massive enough to host 
QSOs, since Ida Angela et al.l l|2008l ) for instance estimated 
a QSO host halo mass of around 3 x 10^^ Mq independent of 
redshift and luminosity from the 2dF-SDSS survey. However 
to establish any dependency of the proximity effect signal 
with the host's mass, they are included for the purpose of 
an extreme low mass range. 

In order to test the influence of the large scale environ- 
ment around QSO host halos we additionally construct a 



sample of lines of sight for each redshift where the origins 
are randomly selected. This sample provides a null hypoth- 
esis since any influence of large scale density fluctuations 
averages out if the origins of the lines of sight are randomly 
distributed in the box. 



3.3 Measuring the proximity effect in spectra 

The proximity effect signature is measured in Lya forest 
spectra by first constructing an appropriate cj-scale for the 
observed QSO. In this study we assume the QSO to have a 

and 



Lyman limit luminosity of L^, 



10''" erg Hz"' s" 



the photoionisation rates of our model UVB field are used. 
We have shown in PI that for higher QSO luminosities the 
proximity effect is more pronounced and less affected by 
the density distribution. Therefore we use a low luminos- 
ity QSO to obtain upper limits of the environmental bias. 
Given the oj-scale, the transmission spect ra is then evalu- 
ated for each line of sight in bins of A log uj iDall'Aglio et al.] 
i2008b 'l. Using the mean transmission per bin, the effective 
optical depth in the bin Teff,Qso(A log cj) is calculated and 
normalised to the effective optical depth in the Lya forest 
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unaffected by the QSO's radiation Ten^Lya- The normalised 
optical depth ^ is thus 

g(Aloga;)="-"-'^^°(^^°g'^). (3) 

The imprint of the proximity effect onto the normalised op- 
tical depth for a given oj-scale becomes 

au) = {l + u;{r)f-^ (4) 

ijLiske fc Willigeill200ll ) where /? is the slope of the Lya ab- 
sorber's column de nsity distribution . Throughout this work 
we assume 13 = 1.5 HimelaLlllool). The UVB photoionisa- 
tion rate Fuvb can then be determined using the proximity 
effect strength parameter a 

This parametrisation was introduced by ijDall'Aglio at al.] 
[2 008b , a) in analyses of observed spectra with an assumed 
UVB as reference. Values of a > 1 or a < 1 indicate a 
weaker or stronger proximity effect than the model, respec- 
tively. The measured photoionisation rate of the UVB is then 
determined using the reference value FuvB.rof multiplied by 
the strength parameter a. In our case the strength parameter 
a indicates any deviation of the measured UVB photoioni- 
sation rate from the input value. The strength parameter is 
determined by fitting Eq. [5]to the binned normalised optical 
depth ^(Alogw). In order to exclude any direct impact of 
the host halo on the strength parameter fit, only data with 
logo; < 2 have been used. 

In PI we found the observed ^ to fiuctuate strongly 
around the analytical proximity effect profile. Any fit of Eq. 
[S]is thus biased by these large fiuctuations which arise from 
the presence of strong absorbers along the line of sight. 

In order to obtain a smoother ^-profile, the wavelength 
scale of the spectrum is transformed into the log u-scale and 
a boxcar smoothing (also known as moving average) with 
the size of A log uj is applied to the transmission spectrum. 
Analogously to the method given above, we then calculate 
the effective optical depth in each smoothed pixel and deter- 
mine the normalised optical depth This results in smooth 
^-profiles allowing the identification of areas dominated by 
strong absorption systems. 

In Fig. [6] two examples of mean ^-profiles taken out 
of the 20 most massive halos sample at z = 4 are given. 
A A log a; = 0.5 is found to yield a good balance between 
smoothing and retaining structure in the ^-profile. For the 
mean ^-profiles shown in Fig. [S] 100 lines of sight have been 
constructed randomly around the DM halos. The two exam- 
ples illustrate a halo without any signs of intervening strong 
absorption systems (left panel) and one where a strong den- 
sity feature is located near the host halo at logt*; — 1.5 (right 
panel). From the halo without strong nearby systems it be- 
comes evident that the profile closely follows the analytical 
one. A fit to this smooth ^-profile is now very robust, since 
the functional form of the profile is well constrained and the 
fit of Eq. [S] is not solely determined by a small number of 
data points. The same is true for cases with strong interven- 
ing absorption systems, as long as the underlying smooth 
profile of Eq. 0] is still visible. Strong deviations from the 
analytical form however, as presented in our second exam- 
ple, are large enough to prevent a clear identification of the 



smooth analytical profile and the obtained strength param- 
eters have to be treated with caution. 



4 RESULTS 

4.1 Null hypothesis: Random locations 

In order to assess whether large scale over-densities affect 
the proximity effect profile and the associated strength pa- 
rameter, we first establish results that are unaffected by such 
large scale density features. A sample of 500 randomly se- 
lected lines of sight originating at random points in the sim- 
ulation will serve as a null hypothesis. No noise is added 
to the spectra. However we will use this sample in the next 
section to discuss the effect of detector noise on our results. 

In Fig. [7| we discuss results obtained for redshifts z = 
4.8, and 3. For each line of sight, the normalised optical 
depth 5 was calculated using the method described in Sect. 
13.31 For each redshift we determine the mean and median ^ 
profile, as well as the ^ probability distribution as a function 
of iu. The ^ probability distribution is shown in colour coding 
in Fig. El 

At z = 4.8 the mean profile follows the analytic prox- 
imity effect model very well, with just a slight increase in its 
slope towards higher lj values. Note however that at large 
ij values of logcj > 2, ^ can only be poorly determined due 
to a very small number of pixels contributing to the Auj 
bins. The median profile as well follows the input model up 
to logo; ~ 0.5. However at logo; > 0.5 the median profile 
steepens strongly and starts to deviate from the input model 
and the mean profile. This indicates a growing asymmetry in 
the ^ distribution with increasing uj (approaching the QSO), 
resulting in the growing discrepancy between the mean and 
median profiles. Considering the ^ distribution function, the 
increasing skewness of the distribution becomes apparent. 
For logu < 1, the width of the distribution stays constant 
and appears symmetrical and normally distributed in the 
logarithmic scale. This indicates a log normal distribution 
of ^ values. However at larger uj values the distribution starts 
to widen up and the peak in probability shifts towards lower 
values of log ^, moving away from the expectations of the an- 
alytical formalism. 

For redshifts z — i and 2: = 3 a similar picture emerges. 
The 5 distribution widens as a function of redshift, with an 
increase in its variance with decreasing redshift. The mean 
profile however always regains the input model well. How- 
ever the discrepancy between the median profile and the 
mean at high uj values increases with decreasing redshift. 
The median profile becomes steeper and steeper, indicating 
a growing skewness of the ^ distribution at high uj. However 
at low cj, the ^ distribution stays symmetric and continues 
resembling a log normal distribution up to logo; ~ 0. At 
z — 2, the ^-distribution shows a strongly increasing vari- 
ance, dominating over the signal of the input model com- 
pletely. We will therefore not consider results from z = 2 in 
this work. 

The increase in variance is dominated by two factors. 
On the one hand the universe evolves and the growth of 
structure increases with decreasing redshift. This introduces 
stronger density contrasts between under- and overdense re- 
gions. On the other hand, the IGM becomes more trans- 
parent with decreasing redshift due to the increasing UV 
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Figure 6. Mean normalised optical depth profiles obtained using a moving average smoothing with a kernel size of A log(tj) = 0.5 for 
two halos taken from the 20 most massive halo sample at redshift 2 = 4 (black solid line with grey shading). The mean halo profiles are 
determined using 100 lines of sight and the shaded area gives the Icr standard deviation. The analytical proximity effect model is given 
as the smooth black line, while the fitted profile is indicated by the dashed line. The left panel shows a halo without strong intervening 
density features and closely following the analytic model. The right panel illustrates a halo having strong density features in its vicinity, 
causing strong deviations from the analytical form. 



background and cosmic expansion reducing the mean den- 
sity of the universe. Therefore at low redshifts, the Lya for- 
est traces denser structures than at high redshift. 

From the simple test of this subsection we conclude that 
with randomly selected lines of sight of random origins, the 
input proximity effect model can be regained with the mean 
^ profile. In our sample this is valid for z ^ 3. The median 
profile however deviates more and more from the mean pro- 
file with decreasing redshift and cannot be used to measure 
the UV background from the proximity effect with the cur- 
rent analytical model, which is formulated for a mean IGM. 

4.2 Dependence on the spectra signal to noise 

We now want to discuss the effect of detector noise on the 
mean proximity effect profiles. Using the 500 lines of sight 
from the null hypothesis sample, noise is added to each spec- 
trum using a signal to noise (S/N) of fO, 20, 50, 100, and 
150. For each S/N sample we determine the mean ^ profile 
and normalise it to the S/N= 150 results. The infiuence of 
noise on the proximity effect signature is shown in Fig. [HI 
where the normalised ^ profiles are shown for z = 4.8 and 
a = 3 as a function of the signal to noise. 

The largest influence noise has on the proflle is at high 
oj values where only a small number of pixels contribute to ^. 
However for uj values below logoj ~ 1.5, the strong infiuence 
of the noise on the profile disappears. At z = 4.8 an increase 
of noise effects can be noted at logoj < 0. There the results 
with a S/N of 10 produce deviations of up to 5% from high 
signal to noise spectra. However already with a signal to 
noise of 50, the ^ profile is regained with an accuracy of less 
than 1%. At 2 = 3 convergence with the low noise results for 
logo; < 1.5 is already achieved with a signal to noise of 20. 
In order to resolve the profile within 1% accuracy at higher 
UJ values, a signal to noise of 50 is needed. 

One has to keep in mind though that these results only 
apply to the combined sample and not to individual lines 
of sight. However this shows that with our approach of not 



considering noise in the spectra, high signal to noise results 
are very well approximated. 

4.3 Large scale environment 

We now want to characterise the large scale density and ve- 
locityenvir onment arou nd our QSO host halos. It has been 
shown by iPrada et al.i (i2006h and iFaucher-Giguere et al.l 
(j2008l ) that the mean density profile around massive halos 
does not reach the mean cosmic density before a halocen- 
tric distance of comoving Vc ~ 10 Mpc. For smaller 
radii the mean DM density around halos can be up to 
10 times above the mean density for halocentric distances 
of comoving Vc ~ 1 Mpc. For even smaller radii the 
density profile is governed by the density profile of the 
host halo. A similar behavio ur has been observed by PI. 
iFaucher-Giguere et al. I ([2008') have further determined the 
bias caused by Doppler shifting of absorption lines due to 
infalling gas in proximity effect measurements. The addi- 
tional bias is found to be similarly important to the effect 
of overdensities. This effect is naturally included in our syn- 
thetic spectra taken from the simulations. 

These large scale features are caused by large scale 
modes perturbing the density distribution at scales of sev- 
eral Mpc. These large modes still experience linear growth 
and increase with decreasing redshifts. Since halos with very 
large masses form preferably where these large modes show 
a positive amplitude, an increase in density above the mean 
cosmic density is expected around very massive halos. Since 
low mass halos also form in regions where there is less power 
on large scales, the density profiles are expected to be less 
infiuenced. 

In Fig. [9] the smoothed and radially averaged overden- 
sity profiles of our halo sample are shown as a function of 
distance and halo mass. A consistent picture to previous 
findings emerges. At 2 = 4.8 the profiles show a steep decline 
in density up to comoving distances of rc ~ 0.4 Mpc 
for the massive halos. This region directly shows the den- 
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Figure 7. Frequency of the normalised optical depth ^ using A logo; = 0.5 calculated from 500 random lines of sight drawn from our 
simulation box, as a function of tJ. The origins of these lines of sight were chosen not to be centred on any specific haloes, but on random 
points in the box. The left panel gives results for z = 4.8 and the right panel 2; = 3. The dashed line marks the input model used for 
generating the spectra including the proximity effect. The black solid line marks the mean ^-profile, whereas the dotted black line marks 
the median profile. 
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Figure 8. The influence of the signal to noise ratio on the mean normalised optical depth ^ for the null hypothesis. Shown are results at 
redshifts z = 4.8 (left panel) and 2; = 3 (right panel). The g-proflles have been normalised to results using a signal to noise of 150. The 
dotted line denotes a S/N of 10, the dashed line a S/N of 20, the dash dotted line a S/N of 50, and the dash triple dotted line a S/N of 
100. 



sity profile of the host halo. For larger distances the slope 
with which the density decreases becomes shallower and a 
transition between the host halo profile and its surrounding 
density environment is seen. The profile then reaches the 
mean cosmic density at a distance above Tc > 10 Mpc, 
independent of halo mass. 

For 2 = 4 the picture is similar, however differences 
between the low mass halos and the more massive ones in- 
crease. The most massive halos themselves now leave their 
imprint up to a radius of rc ~ 1 Mpc, whereas the 
direct influence of the low mass host halos ends at Vc ~ 
0.3 /i"^ Mpc. The low mass halos show a less pronounced 
large scale overdensity with a shallower decline in density 
than the more massive halos. The transition point where 



the density proflle reaches the mean cosmic density is again 
found at a radius of around rc ~ 10 h'~^ Mpc. Furthermore 
the large scale density around the halos shows slightly larger 
values than a.t z — 4.8, which is due to the ongoing structure 
formation and the linear growth of the corresponding large 
wavelength modes. 

At 2: = 3 the differences between the various halo masses 
further increases. The transition point where the profile 
reaches the mean cosmic density starts to show a depen- 
dency on halo mass. The 20 most massive halos reach the 
mean density at a radius of rc > 20 Mpc, whereas the 
low mass halos already reach it at rc > 7 Mpc. Further- 
more the large scale overdensity has grown in density. We 
therefore expect any infiuence of this large scale overden- 
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Figure 9. Mean density environment around all the hales in each 
mass bin. The density profiles are smoothed for better visibility. 
Further the corresponding u) scale is given as reference. The black 
solid line denotes the mean profile around the 20 most massive 
halos, the blue dash-dotted line denotes the mean around the 
A^halo = lO^'^ MfTi halos, and the red dashed line marks the mean 



10^^ Mq halos. The standard deviation is 
= 10^^ Mq halos exemplary for the other 



'-'0 ■ 

around the Mij^io = 
given for the Mhaio 
environments. The first panel shows z = 4.8, the second panel 
2 = 4, and the last panel 2 = 3. The vertical dotted line marks 
the area of log u) > 2 which was excluded in the determination of 
the strength parameter. 
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Figure 10. Mean infall velocity around all the halos in each mass 
bin. The infall velocity profiles are smoothed for better visibility. 
Further the corresponding u) scale is given as reference. The black 
solid line denotes the mean infall velocity around the 20 most mas- 
sive halos, the blue dash-dotted line denotes the mean around the 
Afhalo = lO^'^ Mq halos, and the red dashed line marks the mean 
around the Mhaio = 10^^ ^0 halos. The standard deviation is 
given for the Mhaio = lO-'^'^ Mq halos exemplary for the other 
environments. The first panel shows z = 4.8, the second panel 
2 = 4, and the last panel 2 = 3. The vertical dotted line marks 
the area of log oj > 2 which was excluded in the determination of 
the strength parameter. 
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Figure 11. The a parameter distribution as a function of halo 
mass and redshift. Each black dot represents the a parameter 
obtained using the mean halo profile of 100 lines of sight. For 
each redshift the first line shows results from the 20 most massive 
halos sample, the second line represents the Mhalo = 10^^ Mq 
sample, and the third line the Mhaio = 10^^ Mq sample. The 
bold black lines indicate the sample's mean a parameter, while 
the grey shaded areas indicate the corresponding Icr standard de- 
viation. The solid line marks the fiducial model, while the dashed 
lines indicate deviations of a factor of two in the UVB measure- 
ments. 



sity on the proximity effect measurement to increase with 
decreasing redshift. Additionally we do not expect a depen- 
dency of a possible bias with halo mass at high redshifts, 
since no large differences between the density profiles are 
seen. However at low redshifts, differences due to stronger 
large scale overdensities with increasing halo mass may leave 
an imprint in the proximity effect profile. 

For completeness we show the mean infalling velocities 
as a function of distance to the host for each halo mass bin 
in Fig. [TOl Analogously to the density environment, the in- 
fall velocities grow with decreasing redshift due to the onset 
of structure formation. The infall velocity profiles for the 20 
most massive halos and the A/haio = 10^^ Mq sample be- 
have similarly at all redshifts. Moving away from the halo, 
the infall velocity rises until it reaches its maximum at the 
position where the density profile shows a transition between 
the host halo and its surrounding density environment. At 
larger radii the infall velocity steadily decreases until no sys- 
tematic infall towards the halo is seen anymore. The radii 
where the systematic infall vanishes is about the size of the 
large scale overdensity itself. In the case of the low mass halo 
though, the infall velocity remains constant at low velocities 
over the whole large scale overdensity. 



4.4 Proximity effect strength as function of halo 
mass 

We now want to test whether the proximity effect strength 
parameter a shows any dependency with the host halo mass. 
For each halo the mean ^ profile is calculated from 100 lines 
of sight centred on the halo position. Then Eq.[S]is fitted to 
the mean profile and a proximity effect strength parameter 
a is obtained for each halo in our sample. Any strength pa- 



rameter with log a > results in an overestimation of the 
UV background photoionisation rate in real measurements, 
whereas a log a < resembles an underestimation. 

In Fig. [TT] each halo's mean strength parameter a is 
marked as a function of redshift and halo mass. The a pa- 
rameter distribution is thus obtained for each halo mass bin, 
and its mean and standard deviation are marked by the ver- 
tical bold line and grey area in the plot. At each redshift, 
the mean values of the distribution show no clear depen- 
dency with the host halo mass within the la fiuctuations. 
For redshift z = 4.8, the different halo bins yield similar 
mean values and standard deviations. Considering the vari- 
ance in the distribution, the input model is regained for all 
the halo mass bins at z = 4.8. 

This picture does not change with redshift. Against the 
expectations motivated by the overdensity profiles discussed 
in Sect. 14.31 the data shows no dependency on the halo mass. 
Even at 2 = 3, where the mean large scale overdensity was 
found to be lower around lower mass halos than around high 
mass ones, no significant dependency is seen. The mean val- 
ues are thus consistent with the input model within la. 
The halo to halo variance increases with decreasing red- 
shift and the strength parameter distribution broadens. At 
z = 4.8 the standard deviation of the whole halo sample is 
(j(loga) — 0.08, while at z = 4 we find cr(loga) = 0.20 and 
(j(loga) = 0.36 for 2 = 3. A si milar increase in the varia nce 
has been previously observed bv lDall'Aglio et all Hooii) for 
individual lines of sight of a high-resolution QSO sample. 

From these results we conclude that the proximity effect 
measurements do not show any host halo mass dependent 
bias. The variance in the signal between different halos in- 
creases with decreasing redshift and tighter constraints are 
obtained with higher redshift objects. However even at low 
redshifts, the mean strength parameter regained the input 
value within a factor of less than two. 



4.5 Influence of the large scale overdensity 

In the previous section we have shown that the proximity 
effect strength parameter is not affected by the host halo 
mass. However a large scatter around the mean strength 
parameter is seen in each mass bin, which increases with 
decreasing redshift. We now want to establish the physical 
cause of the scatter, and therefore check if this scatter is 
connected with the density enhancement seen on large scales 
of up to comoving radii of rc = 10 /i~^Mpc. 

For each halo, the mean overdensity profile is calculated 
using the density distribution along each line of sight. Then 
the mean over-density (5)^ between the radii ri and r2 is 
calculated from the mean density distribution around the 
halo using 



1 



r2 — ri 



S{r) dr. 



(6) 



We choose ri = 1 Mpc comoving in order to exclude the 
infiuence of the host halo on the density profile, since we are 
only interested in the large scale density distribution. The 
integral extends to comoving r2 — 10 Mpc which is the 
characteristic size of the large scale density enhancements 
(see Sect. 113]). 

In Fig. [12] the strength parameter a is correlated with 
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Figure 12. Proximity effect strength parameter as a function of 
the mean overdensity in a shell of 1 to 10 Mpc~^ comov- 
ing radius around the 20 most massive halos (black asterisks), 
jW^halo = 10^^ Mq halos (blue triangles), and Mhaio = lO-"^-*^ Mq 
halos (red diamonds). The black dashed line marks a = 1. The 
solid line represent linear regressions obtained from the dataset 
where the dotted lines mark the Ict upper and lower boundaries. 



the mean overdensity in the comoving radial interval = 
[1, 10] h'^ Mpc. The various symbols indicate the different 
halo mass samples. At all redshifts, a correlation between the 
mean density in the vicinity of the QSO and the strength 
parameter of the mean ^-profile is seen. A linear regression 
logo = log (1 + (5)g + a log (1 + 5)ioMpc derived from the 
data, where log (1 + <5)q is the normalisation point and a is 
the slope of the regression. The resulting fit parameters are 
listed in Table[3]together with the Pearson product-moment 
correlation coefficient r of the data set. 

Looking at the z = 4.8 results, the data points align 



Table 3. Table of linear regression results of the overdensity- 
strength parameter relation, together with the data's Pearson 
product-moment correlation coefficients r. The uncertainties in 
the fitting parameters have been determined with the bootstrap 
technique. 



z 


log(l + 5>0 


a 


"^Pearson 


4.8 


-0.15 ±0.07 


0.76 ± 0.08 


0.74 


4.0 


-0.19 ±0.13 


1.07 ±0.16 


0.68 


3.0 


-0.28 ±0.19 


1.11 ± 0.17 


0.60 



clearly on a power law relation. With a Pearson correlation 
coefficient of rpoarson = 0.74 the data exhibits a tight cor- 
relation between the two quantities. The higher the mean 
density around a halo becomes, the stronger the proxim- 
ity effect is biased towards larger strength parameters. At 
z = 4.8 a strength parameter of a = 1 (resembling an unbi- 
ased measurement of the photoionisation rate) is obtained 
for a mean overdensity of (1 + 5)iQMp;. ~ 1.6 ± 0.35. It is in- 
teresting to note that there is no strong segregation between 
the various halo mass bins along the density axis. However, 
the lower mass halos lie to slightly lower overdensities than 
the more massive ones. All the different mass bins cover a 
similar range of large scale densities. Again this indicates 
that there is no distinct connection between the halo mass 
and the mean overdensity on these large scales. Hence, only 
the amount of matter in the greater vicinity of a halo is re- 
sponsible for the large scatter in the strength parameter a 
seen in the previous section. However, a small scatter in a 
of ~ 20% around the linear regression remains at 2 = 4.8. 

Similar findings apply to the lower redshift results. The 
a parameters are still related through a power law with the 
large scale overdensity, however the strength of the corre- 
lation drops to rpcarson = 0.60 at 2 = 3. Comparing the 
z = 4.8 results with the lower redshift ones shows a slight 
steepening of the relation from a = 0.76 ±0.08 a.t z — 4.8 to 
a = 1.11 ± 0.17 at z = 3. The scatter around the power law 
increases to « 35% for z = 4 and 55% for z = 3. Further 
the scatter seems to spread with increasing density, how- 
ever a larger halo sample from a larger simulation would be 
needed to conclusively determine this increase in the vari- 
ance. Again no bias in the strength parameter is found for 
an overdensity of (l + '5)ioMpc ^-^ ^ ^.52 at z — 4 and 
(l + ^)ioMpc~l-8±0.77at2 = 3. 

These results show that the large scale distribution 
of matter around a host halo affects the proximity effect 
strength parameter and biases the measurements of the UV 
background photoionisation rate. There is a power law cor- 
relation between the mean density surrounding a QSO host 
and the a parameter. The bias on the UV background pho- 
toionisation rate can be determined, if the mean overden- 
sity around a QSO is inferred. However, density measure- 
ments from the Lya forest are degenerated with the UV 
background photoionisation rate. The properties of the UV 
background have to be known in order to convert optical 
depths into densities. 

Observationally, iD'Odorico et al.l (|2008l ) have shown 
that QSOs at z = 2.6 are situated in regions showing a 
density excess on scales of 4 proper Mpc. Due to the depen- 
dency of the results on the UV background photoionisation 
rate, their results only constrain the mean density of the 
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large scale overdensity to a factor of about 3. Nevertheless 
we estimate a mean density in a radius of 4 proper Mpc from 
their results. We obtain an upper value on the mean over- 
density of (1 + 5) ~ 3.5 and a lower value of (1 + 5) ~ 1.25. 
According to our results at z = 3, this would translate 
into a bias of the strength parameter of a = 2.1 for an 
overdensity of 3.5 and a = 0.7 for the lower value on the 
mean density around a QSO. Remember though that our 
results only apply to a QSO with a Lyman limit luminos- 
ity of I/i^LL ~ 10"^^ ^'^g Hz~^ s~^. According to PI the bias 
becomes smaller for higher luminosity QSOs. 



5 CONCLUSION 

We have studied the effect of the host halo mass and the 
large scale density distribution on measurements of the prox- 
imity effect strength parameter. The strength parameter 
is observationally used to infer the UV background pho- 
toionisation rate. From a 64 Mpc sized dark- matter 
simulation, we picked the 20 most massive halos, 50 ha- 
los with a mass around A/haio = 10^^ ^^Qj '^nd another 
50 halos with a mass around Afhaio ~ 10^^ Mq at red- 
shifts z = 4.8, 4, and 3. Each halo is individually assumed 
to host a QSO with a Lyman limit luminosity of Lj^^l ~ 
10^^ erg Hz"^ s'^ 

Around each halo, 100 random lines of sight of the den- 
sity and velocity distribution are obtained. Using models of 
the IGM which we calibrated to observational constraints, 
the neutral hydrogen fractions and gas temperatures along 
each line of sight are inferred. To include the increase in 
ionisation of an additional QSO radiation field, the neutral 
hydrogen fractions are decreas ed proportional t o geometric 
dilution of the QSO flux field (|Partl et al.ll2010l ). Then, for 
each line of sight, a Lya forest spectrum is computed. From 
the spectrum we derive the mean proximity effect signa- 
ture with observational methods and measure the proxim- 
ity effect strength parameter a. Since the UV background 
photoionisation rate used in generating the Lya spectra is 
known, the strength parameter quantifies the over- or un- 
derestimation of the UV background. 

In order to assess whether the QSO host environment 
affects the proximity effect signal, a null hypothesis test was 
performed. From the simulation box, 500 lines of sight orig- 
inating at random positions with random directions have 
been obtained. On each line of sight, the proximity effect 
was modelled. From these lines of sight, we obtained the 
distribution of normalised optical depths ^ as a function of 
normalised distance to the QSO ui. From the ^-distribution, 
the mean and median proximity effect signals were deter- 
mined. The mean profiles match the analytic model very 
well and the model UV photoionisation rates are regained. 
At large halocentric distances, the median and the mean pro- 
files match. However the median increasingly deviates from 
the analytic model and the mean profile with decreasing 
distance to the QSO. This indicates that the ^-distribution 
closely resembles a log normal distribution at large distances 
from the QSO, however it is increasingly skewed when near- 
ing the source. Not only is the variance of the ^-distribution 
found to increase with decreasing redshift, but its skewness 
increases as well. 

We further determined the mean density distribution 



and infall velocity around all the halos in each mass bin 
as a function of redshift. For all redshifts the density pro- 
file does not immediately reach the mean cosmic density 
at the border of the dark matter halo, however it stays 
above the mean density up to comoving halocentric radii 
of > 10 Mpc. For redshifts z = 4.8 and 2 = 4 the 
size of this large scale overdensity is independent of the halo 
mass. However higher mass halos show in the mean a larger 
overdensity than lower mass halos. At redshift z = 3 the size 
of the large scale overdensity is smaller for the lower mass 
sample than for the more massive ones. The velocity field 
shows a mean infall up to distances of 30 Mpc around 
the halos, reaching for the most massive halos 80, 100, and 
150 km/s for redshifts z = 4.8, 4, and 3, respectively. Over- 
density and infall velocities act together in the derived bias 
of the proximity effect. 

The mean strength parameter per halo mass bin does 
not show a dependency with the halo mass. For each halo 
mass bin, a mean strength parameter which is consistent 
with the input model is regained within la standard devia- 
tion. However the halo to halo variance is found to increase 
with decreasing redshift from cr(loga) = 0.08 at z — 4.8 to 
cr(loga) = 0.36 at z = 3. 

The strength parameter is found to correlate with 
the mean density measured in a shell of comoving 1 — 
10 h~^Mpc. We fit a power law to this correlation. Regions 
with a mean density below the cosmic mean show a stronger 
proximity effect than regions with densities above the cos- 
mic mean. The correlation is tightest at z = 4.8 and the 
scatter around the power law increases with decreasing red- 
shift. Furthermore, the power law is found to steepen slightly 
with decreasing redshift. From these results we find that an 
unbiased UV background photoionisation rate can only be 
obtained if the mean density around the QSO is between 
(1 + ^)ioMpc ~ 1-5 ± 0.52 and 1.8 ± 0.77. 

If a possibility exists to determine the mean density 
around a QSO, the bias arising from the large scale over- 
density can be corrected for. However density measurements 
from Lya forest spectra are degenerate with the UV back- 
ground photoionisation rate. Due to this degeneracy, previ- 
ous measurem ents of the density distribution around QSOs, 
such as D'Od orico et al.l (2008), have only been able to con- 
strain the mean density up to a factor of 3. By inferring 
the bias of the proximity effect at z ~ 3 according to their 
density measurements, the strength parameter ranges from 
a — 0.7 to a = 2.1 for a QSO with a Lyman limit luminosity 
of L^LL = 10^' erg Hz'^ s'V 
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